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Abstract 

On the base of symplectic quantum tomogram we define a proba¬ 
bility distribution on the plane . The dual map transfers all observ¬ 
ables which are polynomials of the position and momentum operators 
to the set of polynomials of two variables. In this representation the 
average values of observables can be calculated by means of integration 
over all the plane. 


1 Introduction 

In the twentieth century there were several attempts to describe quantum 
mechanics not with help of operators, but by the language of functions, as 
in classical mechanics. For this it is necessary to assign a certain (maybe 
quasi) probability distribution function to the quantum state, and another 
function dehned on the same space to the operator of observable. In other 
words, let p and a be a density operator (a positive Hermitian operator with 
unit trace) and an observable, respectvely. Then, these operators should be 
associated with the function p{X) and a{X), 

p tA p(A), a, tA ci(X), X G (<T, p), (1) 
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on the measurable space (A”, p), wherein the average value of a in the state 
p can be calculated by the formula 


<a>p= p{X)a{X)dp{X). (2) 

Jx 

The realisation of this program resulted in the requirement of taking the 
appropriate set of test functions {p(X)} for which {a(A)} are distributions 
(generalised functions). 

In 1932, Wigner [1] introduced the function, later named in his honor 


W{q,p) 


1 


p(9+ |)e 


where p(-,-) is the density matrix of p in the coordinate representation. 
Wigner noticed that for functions either only of the position operator q or of 
the momentum operator p alternatively, the average values can be calculated 
according to the formulas 


< F{q) >~p= / F{q)W{q,p)dqdp 
Jr'^ 


<G{p)>p= G{p)W{q,p)dqdp 

JR2 

In 1949 Moyal [2] generalized this result to the case of arbitrary poly¬ 
nomials of q and p. He showed that the mathematical expectation of the 
symmetrized polynomial of the form 

n 

= (3) 

can be calculated by the formula 

< {q'^p"'}s >p= f Fmn{q,p)W{q,p)dqdp 


at that Fmn{q,p) = q^p'^- Thereby Moyal proposed method of calculating 
average values for a broad class of observables through the Wigner function 
along the lines of classical mechanics. However, it should be noted that the 
Wigner function W{q,p) is not positive in general. So it can not be regarded 
as a classical probability distribution function. 

In [31llj the optical quantum tomogram w(X, p) was introduced, dehned 
as the Radon transform of the Wigner function 

w(W, (/?)= / W{q,p)6{X — qcos p — psmp)dqdp. (4) 

iR2 
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Function w(X, ip) is non-negative everywhere and a one-parameter set 

{w(X,9?), ip e [0,27r)} 

consists of a probability distributions on the real axis. 

In [S] it was noticed that the optical tomogram w^(X, ip) of a quantum 
state p can be calculated directly from the density operator 

Wp(X, p) = Tr(p(5(XJ — q cos p — psin p)). 

Similarly there was introduced a more general notion of symplectic quantum 
tomogram 

ojp{X,p,p) = i:i{p5{Xi - pq-vp)). (5) 

Symplectic tomogram has the property of homogeneity 

a;(AX, A/i, Xu) = —a;(X, p, v) 

and it contains the same information as a density operator. The tomographic 
representation of quantum mechanics was developed in many branches during 
elapsed time Such new approaches as C'*-algebras formalism appeared 

to be connected with this topic [8]. 

In piTn] it was studied the following question. Suppose that one has taken 
for the correspondence p p{X) in ([1]) the map from the set of states with 
sufficiently smooth density matrices p(-,-) in the coordinate representation 
to the set of optical tomograms w{X,p) (in the case, variable X = {X,p)). 
Then, the dual map a —?■ a{X) in ([T]) is given by the extension of 

a{X, (p) = —2 tt \im.Tr{a{X + ie — cos{(f))q — sin(^)p)“^) (7) 

£—>■0 

from the set of density operators a. It was shown that for the symmetrized 
product of canonical quantum observables dehned in ([3]) the dual 

map results in 

amn{X,p) = X^+^H^^^{p) ( 8 ) 

where functions are biorthogonal to the system {cos"'”''™'”^ psin^ p}, 

i. e. 

2n 

j cos”+™"^((p) Sin^((p) Hn+m-s,s{‘f)dy^ = ^ks- (9) 

0 

In the present paper we pass from the polar coordinates {X, p) to the 
Cartesian coordinates {x,y). Our aim is to introduce a probability distribu¬ 
tion 0(a:, y) on the plane such that the average values can be calculated 
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by the formula 


+ 00 +00 

<{q'^p'^},>p= j j fi,{x,y)VL{x,y)dxdy (10) 

— 00 —00 

where function f^{x,y) corresponds to in the Cartesian coordi¬ 

nates. 

2 The comparison with the traditional ap¬ 
proach 

To calculate the average value of the observable d in the state p in the 
traditional representation of quantum mechanics one should take a trace m 

< a >p= Tr{pd). (11) 

In this section we discuss the advantages of applying (|2]) instead of the tra¬ 
ditional approach (ITT]) . 

First and foremost, the optical tomogram (jl|) is a measurement of the 
homodyne quadrature [121 [13] • After one obtains it, the density operator 
p should be reconstructed from the experimental data. If we consider the 
tomogram as the main unit describing the quantum system, instead of the 
density operator, then this reconstruction is not needed. It allows one to 
take maximum advantage of experimental data. Such an approach results in 
a number of signihcant tasks concerning the situation in which one should 
reconstruct a tomogram in the case of incomplete information [I3HIZ|. 

On the other hand, the tomographical technique can be successfully used 
for comparing classical and quantum systems in the sense of [I8]. Indeed, 
the role of a density operator in the classical statistical mechanics is played 
by a probability distribution on the phase plane. Therefore, to compare 
classical and quantum cases probability distributions are preferable to density 
operators. In this framework, one can discriminate what kind of probability 
distributions on the plane refer either to the classical or the quantum case. It 
is worth mentioning that using the symplectic tomogram allows one to take a 
limit for the centre-of-mass tomogram [12] under the condition that a number 
of particles tends to inhnity [23] • Surely these are but several examples of 
the possibilities of using tomographic picture. 
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3 Tomography on the plane 


In [9l[Tn] the optical tomogram w{X,(j)) (jl]) plays a role of a probability 
distribution on the set M x [0,27r]. Indeed, 


27r +00 


w(X,(p) > 0, — 



w(X, ip)dXd(p = 1. 


( 12 ) 


0 —oo 


In turn, the information containing in the set {w(X, (p)} is redundant due 
to that 

w(—X,ip)=w{X,ip + 7r) (13) 

(it immediately follows from (|4])). It results in the knowledge of {w(X, (j)), X > 
0, 0 G [0,27r]} allowing reconstruction of the entire tomogram. The pair 
{X, (j)) G M+ X [0,27r] can be considered as the polar coordinates on the plane 
Thus, it looks natural to change the polar coordinate to the Cartesian 
coordinate. This should give a transformation of w(X, (p) to a probability dis¬ 
tribution Q{x, y) on the plane Let us dehne a function r2(x, y) in terms of 
the symplectic tomogram (|5]) by the formula (here we used the homogeneity 

O) 

1 / rv;- X y 


Vt{x,y) = 


It follows then that 


^/x"^ + y' 


-.u 


\/x‘^ + 


y 


x"^ + y"^'' x"^ + y"^ 


(14) 


= u{x^ + y‘^,x,y). 


^{x,y)>0, — / il{x,y)dxdy = 1. 


IT 


(15) 


Moreover, let us assume that a{X, cp) is the symbol of observable d in the 
optical tomography representation ([7]). Let us presume that 

fa{X cos{(j)), X sm{(j))) = a{X, (p). 


Assertion 1: 


<a>p= / fa{x,y)n{x,y)dxdy. 


Proof . Consider the right-hand side of flTU]) and pass to the polar coordi¬ 
nates 


+ 00 +00 



fa{x,y)n{x,y)dxdy = 


— OO —oo 
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27r +00 



3 {x^ if] 

fa{r cos if, r sin 99)f2(r cos yj, r sin — -dr dip = 


0 0 

27r +00 


a(r, <p) 



fair cos (/?, r sin ip)-uj{r, cos 9?, sin ip)rdrd(p = 
r 


0 0 


As required. 
Put 


27 r +00 



fair cos ip, r sin y9)w(r, ip)drdip. 


0 0 


f^n{.X cos ip, X sirup) = am„(A:,9?), 


( 16 ) 


where amn is defined in ([8]). The function fmn{x,y) determined by ffThjl is a 
symbol for the observable d = {q"^p‘^}s- 

Let us consider the Hilbert space consisting of all polynomials of A^th 
degree on the plane equipped with the scalar product 


if, 9 ) 




R2 


y) f(x,y)g(x,y)dxdy. 


(17) 


Assertion 2 : The functions ifmnix,y))m+n=N form the biorthogonal sys¬ 
tem to the functions with respect to the scalar product flTT)) . 

Proof . 

Passing to the polar coordinate we obtain that fiX, <h) = X^F{(f)) and 
g{X, $) = X^G{(j)) are biorthogonal iff F{(j)) and G'(0) are bioorthogonal on 
the unit circle. 

In the Appendix we find the functions fmn for m -|- n = 2 and m -|- n = 3. 
Thus, 

f 2 oix,y) = fiSx^ -y^), fuix,y) = -xy, fo 2 ix,y) = fi3y^ - x^), 

2n TT Zn 

2 2 

hoix,y) = -{x^ - xy‘^),f 2 iix,y) = -idx^y-y^), (18) 

TT TT 

2 2 

fi 2 ix, y) = -(5x?/^ - x^), /o 3 (x, y) =-{y^ - x^y). 

TT TT 
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4 Some simple examples 

Let us calculate Q{x,y) determined by flT^ for coherent states. The sym- 
plectic quantum tomogram of a coherent state \a > with the wave function 


< X\a >= 


71 


1/4 


exp 


r- o?' IckP 

2 2 2 


is known to be equal to (see, he., m) 

1 [ {X — \/2Re{a)y — \/2Im{a)u)'^ 


Ua{X, y, u) = - 

Applying ffTTl) we obtain 
1 


: exp 


fJ? + Z/2 




: exp 


+y‘^ — \/2Re{a)x — \/2Im{a)yy 


In particular, for a vacuum state |0 > it implies 

1 


^o{x,y) = 




exp (—— y )■ 


The number of particles operator N = ^ has the symbol 

= ^U2o{x,y) + fo,2{^,y) - ^) = + 

due to flTSl) . It immediately follows that 


< N > 0 = / no{x,y)fj^{x,y)dxdy = 0 


Analogously, for the excited states |r?. > of harmonic oscillator, with the wave 
functions 

< X\n >= -—rr—;^=HniX) exp 

where Hn{X) is the Hermitian polynomial of nth degree, the symplectic 
tomogram equals 



Un{x,yi,u) = 


m 


X 


2"n! ^7r(/x2 + u^) ” \^7i{y^ + 


exp 




y? + z/2 
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Then, it follows from flTT]) that 


Vtn{x,y) = 


2-n! ^7r(a:2 + ?/2 
The same reasoning gives us 


-.Hi 


71 


exp(—— y'^). 


< N >n= / Qn{x,y)fj^{x,y)dxdy = n. 


5 Conclusion 

Based upon the symplectic quantum tomogram we introduced the probabil¬ 
ity distribution on the plane ([HD. Then, we have presented the dual map 
allowing to calculate the average values of quantum observables in the quan¬ 
tum state determined by this distribution. Finally, we have calculated the 
symbols of quantum observables which are polynomials of the position and 
momentum operators of second and third degrees. Few simple examples of 
calculating the probability distribution on the plane and the symbols of ob¬ 
servables are given. We suppose that our approach can serve for solving the 
task of reconstructing a tomogram from incomplete information about the 
system. We shall treat with this problem in the future. 
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Appendix 

Here we find a biorthogonal system f^\x^y) to the auxiliary system 
g^r\x, y) = x^~^y^ with respect to the scalar product 

2 


{f:9) = 


n\ 



e fl''\x,y)g^J^\x,y)dxdy 


Obviously fi^\x,y) = f, f 2 ^\x,y) = f. Look for now fl^\x,y). We 
have: 

9 ?\x, y) = x^, 9?\x, y) = xy, g^^\x, y) = y^ 

Functions f^^\x,y) are polynomials of the second degree relative x and y, 

(‘2'\ 

therefore they represent a linear combination of functions gs ix,y) that is 

fk^\x, y) = ttkx"^ + bkxy + Cky"^ 

While the basis is gi^\x,y) we can write the column coordinate functions 

/ak\ 

fj^\x,y) as j 1 • To calculate the scalar products {f,g) we write the 

\CkJ 

Gram matrix in the basis of gi^\x,y): 

({9?\9?) {9?\9?^) {9?\9?)\ 

r= {g^^\g?) {g?\g?^) 

\i9?\9?^) {9?\9?) {9?\9?) 


where 

(ft'hsf) = 


+00 


+00 



e gf\x,y)g'f\x,y)dxdy = I x^ * ^^dx I y'^^^e dy 


.(2) 


The values of these integrals are known, so 

^3 0 1 ' 


TT 


r = - 0 1 0 
^ u 0 3. 


Condition of biorthogonality for example for/| is written as 

(fdJ?) = 1 
(/h.sf) = 0 
(/f’.sft = 0 
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that is 

( 3ai + Cl = ^ 

I h = 0 

^ oi -|- 3ci = 0 

The solution is y) = — y"^). Similarly, the other elements of the 

system are 

f?\x,y) = -xy, f?{x,y) = ^{-x^ + 3y‘^) 

TT /TT 


For n = 3 we have 

9f\x,y) = gf\x,y) = x^y, gf\x,y) = xy‘^, gf\x,y) = y^ 

Gram matrix in the basis of gi^\x,y) is 

/5 0 1 0\ 
n 0 10 1 
“8 10 10 
\0 1 0 5/ 

By similar calculations we obtain 

f?\x,y) = -{x^ -xy^) 

TT 

f2^\x,y) = -(5x^2/ - y^) 

TT 

= -(5xy2 - x^) 

71 

f4\x,y) = -{y^ -x^y) 

TT 
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